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Abstract 

This paper summarizes the main results of our recent study of the 
non-Blazhko RR Lyrae stars observed with the Kepler space telescope. 
These stars offer the opportunity for studying the stability of the pulsa- 
tions of RR Lyrae stars and for providing a reference against which the 
Blazhko RR Lyrae stars can be compared. Of particular interest is the 
stationarity of the low-dispersion (a < lmmag) light curves constructed 
from ~ 18000 long-cadence (30-min) and (for FN Lyr and AW Dra) the 
~150000 short-cadence (1-min) photometric data points. Fourier-based 
[Fe/H] values and other physical characteristics are also derived. When 
the observed periods are compared with periods computed with the War- 
saw non-linear convective pulsation code better agreement is achieved as- 
suming pulsational L and M. values rather than the (higher) evolutionary 
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L and M values. 

1. Introduction 

In the 105 deg 2 field of the Kepler space telescope there are ~40 RR Lyr stars 
brighter than 17th magnitude being observed at both long and short cadence (i.e., 
integrated exposure times of ~30-min and 1-min, respectively). With the notable 
exception of RR Lyrae itself (Kolenberg et al. 2010, 2011) little was known about 
the stars prior to the launch of telescope. Several of the RR Lyrae stars are c-type 
stars pulsating in the first overtone mode (Moskalik et al. 2011, in preparation). The 
remainder are ab-type stars pulsating in the fundamental mode, approximately half of 
which exhibit the Blazhko effect (see Benko et al. 2010) and show phase and amplitude 
variations of the light curves occurring on time scales of ~10 to hundreds of days. The 
other half of the RRab stars appear to oscillate in a remarkably constant manner and 
are the subject of our recent paper "Fourier analysis of non-Blazhko ab-type RR Lyrae 
stars observed with the Kepler space telescope" (Nemec et al. 2011), the results of 
which are summarized here. Elsewhere in these conference proceedings Kinemuchi 
(2011) reviews the general properties of the Kepler RR Lyrae star observations, and 
Kolenberg (2011) discusses the Blazhko RR Lyrae stars. 

After deriving pulsation periods, times of maximum light, and total amplitudes, 
and characterising the light curves using Fourier decomposition methods, attention is 
turned to [Fe/H] values and how the stars compare with other RR Lyrae stars in our 
Galaxy and in the Large Magellanic Cloud. After establishing that there are systematic 
offsets between the V and Kp Fourier parameters, the Kp versions of the correlations 
established by Kovacs, Jurcsik, Walker, Sandage and others are used to derive physical 
characteristics for the stars. These include the dereddened colour (B — V) , mean 
effective temperature, T e s, surface gravity, log g, and pulsational and evolutionary lu- 
minosities (L) and masses (M). Finally, the results are compared with new theoretical 
model results computed with the Warsaw non-linear hydrodynamics code. 

2. Observations and Data Analysis 

For each of the 19 stars listed in Table 1 the K epler raw-flux data acquired over the 
first 420 days of the telescope's operation (Q0-Q5) was analysed. The mean apparent 
magnitudes of the stars (from the Kepler Input Catalog) range from {Kp) = 12.7 to 
17.4 mag, with distances ~3-23 Kpc in the direction (l,b) = (76.3°, 13.5°). The last 
two columns of Table 1 identify the particular Kepler data that was analyzed. 

Figure 1 shows the raw-flux (left panel) and fully-transformed .fTp-magnitude (right 
panel) as a function of time for one of the non-Blazhko stars, V894 Cyg (graphs such 
as these are typical). The flux is in units of e~/s/cadence (logio scale) and 18911 
long cadence (LC) aperture photometry measurements are plotted. Since the 1.4-m 
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TABLE 1 

Non-Blazhko ab-type RR Lyrae stars in the Kepler field: 
Pulsation Periods, Times of Maximum Light, and Analyzed Data 
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telescope was rolled by 90° every three months resulting in a given star being observed 
with different pixels every 'quarter', it was necessary to correct for sensitivity variations 
from quarter to quarter, and for drift within each quarter. The Kp magnitudes were 
derived from the fully normalized and detrended photometry. The scalloped pattern 
seen at maximum and minimum light, and the moire pattern seen everywhere, result 
from the pulsation period for this star, 0.571 d, not being evenly divisible by the LC 
sampling period, ~30 min, giving rise to 27 observations per pulsation cycle. The 
constancy of the range in the total amplitude is remarkable, as is the continuous moire 
pattern. 

The pulsation period for each star (fundamental mode) was derived from the high- 
precision aperture photometry using the PERIOD04 period-finding program (Lenz & 
Breger 2005) and a version of the CLEAN program written by Dr. Seung-Lee Kim (see 
Nemec, Walker & Jeon 2009). Typically the LC data consisted of ~18000 points, and 
the short cadence (SC) observations of FN Lyr and AW Dra (Q0 and Q5) amounted 
to another ~ 150000 data points per star. The uncertainties in the derived periods are 
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Figure 1. — Raw-flux (left) and i^p-magnitudc (right) long-cadence photometry for V894 Cyg, 
one of the 19 non-Blazhko RR Lyrae stars in the Kepler field. The processing of the raw-flux data 
consisted of normalizing the different (Q0-Q5) sensitivities, removing (minor) trending within each 
quarter, stitching together the quarters, and shifting the magnitudes to (Kp) = 13.293 mag. The 
times are BJD (Barycentric Julian Date) minus 54953.0 (the first day of Kepler observations) . 




Figure 2. — (Left) Light curve for FN Lyr derived from Q1-Q5 long-cadence photometry; (Right) 
Light curve for AW Dra derived from short-cadence Q5 photometry. (Lower) The residual plots in 
both cases are from the 38- term Fourier fits (observed minus fitted). The standard deviations of the 
fits in both cases are below 1 mmag and would be even smaller with more terms in the Fourier fit. 

~l-2xl0 7 d, depending on the particular data set and pulsation period, P. Times of 
maximum light, t , accurate to ±0.0005 d also were computed. The P and t values 
are summarized in columns 4 and 5 of Table 1. 
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Phased light curves for all 19 non-Blazhko stars were plotted using the derived 
pulsation periods and times of maximum light. Two such light curves are illustrated 
in the top panels of Figure 2. The light curve for FN Lyr is based on 18338 LC data 
points acquired over ~400 days in Q1-Q5, and that on the right is for AW Dra and is 
based on 135380 SC data points acquired over 90 days in Q5. When the light curves 
for all the stars are plotted one clearly sees variations in the 'risetimes' (i.e., the time 
from minimum to maximum light, expressed in phase units), the stars with the shortest 
risetimes having the largest amplitudes and a tendency to have a bump on the rise to 
maximum light. 

Also analyzed were new high precision ground-based V photometry for three of 
the stars (NR Lyr, FN Lyr and AW Dra), and All Sky Automated Survey (ASAS) 
V, I photometry (see Pigulski et al. 2009) for nine of the brightest stars. The high- 
precision calculations were instrumental in establishing the V-Kp offsets in the Fourier 
parameters (needed to permit use of the Kovacs, Jurcsik et al. relations); and the 
ASAS-North data (http://www.astrouw.edu.pl/asas/kepler), although not providing 
high-precision Fourier parameters, were sufficient for deriving reliable loop diagrams 
in the HR-diagram, from which it is clear that in every case the non-Blazhko stars are 
bluest when brightest (as expected). 

By combining the historical data (going back more than 100 years for FN Lyr, 
NQ Lyr and AW Dra) with the Kepler data we were also able to derive period change 
rates for some of the stars. In particular, the periods for FN Lyr and AW Dra were 
seen to be increasing. Unfortunately the baseline of the Kepler observations is still too 
short to derive significant period variations directly. 

Finally, one of the non-Blazhko stars, KIC 7021124, was discovered to be doubly- 
periodic and to have properties very similar to V350 Lyr (Benko et al. 2010). Both 
appear to be pulsating in the fundamental and second-overtone modes simultaneously. 
The period ratios are almost identical, P 2 /P = 0.59305 for KIC 7021124, and 0.592 
for V350 Lyr. When masses are derived for these two stars using the Petersen P2/P0 
vs. Pq diagram the best agreement with the models is obtained, in both cases, for a 
high luminosity and a high mass: L/L Q = 70 and M./M.Q = 0.75. 

3. Fourier Analysis 

Fourier decomposition of the light curves was performed on the fully calibrated 
photometry for all the stars by fitting the following Fourier series to the apparent 
magnitudes: 

F 

m(t) =A o + ^Asin[27ri/ o (*-to)+0i], (1) 

i=l 

where m(t) is the apparent magnitude (either Kp for the Kepler data, or my for 
the ground based ^-photometry) , F is the number of fitted terms, fo is the (assumed 
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Figure 3. — Variations in time of the Fourier parameters PHI31, R21, PHI1S and Al derived 
for the typical non-Blazhko star NQ Lyr (left panel) and for the known low-amplitude Blazhko star 
V783 Cyg (right panel) . Each point corresponds to the variable of interest as determined from a single 
pulsation cycle. 



single) pulsation frequency of the star (=P ~ ), * is the observed time of the observation 
(BJD-54953 for the Kepler data), to is the time of maximum light used to phase the 
light curves (so that maximum light occurs at zero phase), and the A; t and <pi are 
the amplitude and phase coefficients for the individual Fourier terms. The Fourier 
calculations were made with two FORTRAN programs, one kindly provided by Dr. 
Geza Kovacs, and the other by Dr. Pawel Moskalik. The assumed pulsation periods 
were those derived either from the Kepler data or from the period change rate analysis, 
and the t values were calculated numerically. 

The Fourier analyses were performed on various subsets of the data. The bottom 
panels of Fig. 2 show residual plots after fitting 38-term Fourier descriptions to all the 
data. In both cases the standard deviations about the mean curve are below 1 mmag, 
specifically ~60 /xmagn, and the greatest differences occur where the slope is steepest. 
For the purposes of using extant Fourier correlations to derive physical characteristics 
only the lower order Fourier terms are needed. Owing to the small risetimes as many 
as 100 or more terms could be needed for a normal distribution of the residuals. 

Another way of looking at the data was to calculate Fourier coefficients and pa- 
rameters for each pulsation cycle. For the long cadence (30 min) data there were ~24 
observations per bin for an RRab star with a period of 12 hours, and fewer (more) 
observations per cycle for stars with shorter (longer) periods. None of the 19 stars ex- 
hibits the recently discovered "period doubling" effect, thus confirming the result found 
earlier by Szabo et al. (2010). For the Fourier calculations 'direct Fourier fitting' (dff) 
rather than 'template Fourier fitting' (tff) methods were used (see Kovacs & Kupi 
2007), and for each star the resulting time series were plotted for four Fourier parame- 



90 



Nemec, Smolec, Benko et al. 



ters: A±, the first term in the Fourier series (see right panel of Fig. 2); <f>i, the phase of 
the first term; R 2 i, the amplitude ratio A 2 /Ai, and the Fourier parameter found 
by Simon, Kovacs and others to be one of the most significant variables for deriving 
physical characteristics. Typical time series are illustrated in Figure 3, on the left for 
the unmodulated RRab star NQ Lyr having intermediate amplitude (Ai=0.279 mag) 
and intermediate brightness, (Kp) = 13.075 mag, and on the right for V783 Cyg, a low- 
amplitude Blazhko star with the shortest Blazhko period (Pb = 27.7 ± 0.04 d) among 
the stars in the sample. For NQ Lyr linear trend lines were fit to each time series and 
in almost every case the slope is zero to within the systematic and random uncertain- 
ties. The mean values for the four parameters are as follows: <ft s zl = 5.0958 ± 0.0023, 
with residual standard error, a = 0.002; R 21 = 0.4710 ± 0.0006, with a = 0.0006; 
(pi = 3.961 ±0.001, with a = 0.001; and A x = 0.28016 ±0.00014 mag, with a = 0.0001. 
Since the other non-Blazhko stars show random variations of the Fourier parameters 
similar to those shown here for NQ Lyr these means and errors provide a measure of 
the typical uncertainties. Inspection of the NQ Lyr time series in Fig. 3 shows that 
the stability over the ~420 d interval (Q0-Q5) is exceptional, not only for each Fourier 
parameter but for the ensemble of parameters. To better illustrate the remarkable sta- 
tionarity of the light curves for all 19 sample stars a set of 'animated gif light curves 
was prepared and these are available in the electronic version of the MNRAS paper. 

In Table 2 some of the resulting Fourier coefficients and parameters are summa- 
rized. Column 2 gives the standard deviation about the mean light curve (see Fig. 2), 
most of which are below one mmag. Column 3 contains the first term in the Fourier 
series, which correlates extremely well with the total amplitude. Columns 4 and 5 con- 
tain amplitude ratios (defined above) and columns 6 and 7 contain phase parameters 
(see Simon 1988 and references therein). 

It is of interest to compare the derived Fourier parameters for the 19 non-Blazhko 
stars in the Kepler field (derived using the Kp photometry but transformed to V 
values using the offsets) with parameters derived for Galactic globular clusters and for 
RR Lyrae stars in the Large Magellanic Cloud. 

Figure 4 compares the non-Blazhko stars with 177 RR Lyr stars in several well- 
studied globular clusters (from Kovacs & Walker 2001). The same symbols and colour 
coding are used in all four panels. The upper-left panel is analogous to the P-A tot 
diagram (see Fig. 6 below) and was used to derive [Fe/H] values for the Kepler stars. 
The two diagonal lines represent mean relations when the globular cluster data (small 
black dots) are sorted into two [Fe/H] bins (the 'BB' subscript refers to the Butler- 
Blanco system, which approximates the Carretta-Gratton system): a metal-poor bin 
consisting of 19 stars (surrounded by blue squares) with metallicities between —1.70 
and —1.99 dex and average —1.8 dex (BB-scale); and an intermediate-metallicity bin 
consisting of 39 stars (circled in red) with [Fe/H]BB between —0.97 and —1.23 dex and 
average —1.1 dex. The equations of the lines are: 0^ = 5.556 logP + 0.2920 (upper 
metal-poor bin) and 0^ = 6.200 logP±3.615 (lower intermediate metallicity bin). Four 
of the Kepler non-Blazhko stars are clearly richer than [Fe/H] B B=— 1-1 dex, while the 
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TABLE 2 

Selected Fourier Parameters and Derived Metallicities for the 
Kepler non-Blazhko ab-type RR Lyrae stars 
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remainder are apparently more metal poor. The upper-right panel of Fig. 4 shows 
that the Kepler stars appear to be drawn from a distribution similar to that of the 
globular clusters. In particular, the Kepler stars with low _R 31 values are not unusual, 
except possibly that they all have relatively high R21 values. Since the globular cluster 
stars of higher metallicity (red open circles) all tend to reside on the right side of the 
diagram this separation is probably a metallicity effect, supporting our conclusion that 
V784 Cyg, V782 Cyg, KIC 6100702 and V2470 Cyg are metal rich. Likewise, the 
lower-left panel shows that the majority of the Kepler stars do not differ from the 
stars in globular clusters. Note too that there is very little metallicity discrimination 
in this plane. The two stars located at the extreme upper edge of the envelope of 
the globular cluster distribution are AW Dra and V784 Cyg. Finally, the lower-right 
panel shows close agreement between the phase parameters of the Kepler and globular 
cluster RR Lyr stars, which supports the conclusion drawn earlier (lower-right panel 
in Fig. 4) that there is a strong approximately-linear correlation between 02i an d 4>li- 
The diagram also shows very little dependence on metallicity. 

Figure 5 compares the R21, R31, 02i an d 03i Fourier parameters for the 19 Kepler 
non-Blazhko RR Lyr stars with those for the field RR Lyr stars in the central regions 
of the LMC. The LMC data are from the massive OGLE-III catalog of variable stars 
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Figure 4. — Four panels comparing the Fourier parameters (transformed to the F-band) for the 
Kepler non-Blazhko stars (large black squares) and for the 177 RR Lyr stars located in several Galactic 
and LMC globular clusters (small black dots). The cluster RR Lyr stars (from Kovacs & Walker 2001) 
are in the globular clusters M2 (12 stars), M4 (4), M5 (12), M55 (4), M68 (5), M92 (5), NGC 1851 
(11), NGC 5466 (6), NGC 6362 (12), NGC 6981 (20), IC 4499 (49), Ruprecht 106 (12), NGC 1466 
(8), Reticulum (8) and NGC 1841 (9). The points for the 19 most metal-poor globular cluster stars 
are surrounded by blue squares, and the points for the 39 most metal-rich globular cluster stars are 
circled in red. 



by Soszynski et al. (2009), which comprises almost 25000 RR Lyr stars. All the 
parameters are V-band values, the LMC values having been transformed from / to 
V using the transformation equations given by Morgan, Simet & Bargenquast (1998), 
and the values for the Kepler RR Lyr stars transformed from Kp to V using the well- 
determined offsets from diagrams such as Fig.6 below. The upper-left panel shows 
that the R21 values for the Kepler stars are typical for ab-type RR Lyr stars and they 
all lie in the narrow range 0.45 to 0.51. The period ranges of the RRab stars also 
are similar. Where the LMC stars differ from the Kepler stars is in the much larger 
R21 range; in particular, none of the Kepler stars are among the LMC stars in the 
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Figure 5. — Four panels comparing Fourier parameters (V-band) for the Kepler non-Blazhko 
ab-type RR Lyr stars (large black squares) and for 24905 RR Lyr stars in the central regions of the 
Large Magellanic Cloud. The LMC data are from the Soszynski et al. (2009) OGLE-III study, with 
colour and symbol coding as follows: red filled dots (17693 RRab stars); blue filled squares (4957 RRc 
stars); green open triangles (986 RRd stars), and purple open diamonds (1269 RRe stars). Because 
of the limited x- and y-ranges not all of the LMC RR Lyr stars are represented in the graphs. 



'dropdown' at logP ~ —0.19. It is suspected that the majority of these are RRab 
Blazhko variables. In the upper-right panel of Fig. 5 the R21 and R31 amplitude 
ratios are compared. This graph puts into perspective the small ranges of _R 2 i an d R31 
seen for the Kepler non-Blazhko RRab stars, both relative to that for the LMC RRab 
stars, and relative to the ranges seen for other types of RR Lyr stars (RRc, RRd and 
RRe). The upper left panel of Fig. 4 has already shown that the Kepler stars sort into 
four metal-rich and 15 metal-poor stars, and that the metal-poor stars are similar to the 
RR Lyr stars found in Galactic globular clusters (both Ool and OoII types); however, 
the metal-rich stars have no counterparts, at least among the Kovacs & Walker (2001) 
sample of globular cluster RR Lyr stars. In the bottom panels of Fig. 5 (both of which 
are metallicity diagnostics) we see that both the metal-rich Kepler RR Lyr stars (the 
clump of four stars at logP ~ —0.3 and §\ x ~ 4.3) and the metal-poor Kepler RR Lyr 
stars (the majority with lower $ c 2l values) have counterparts in the LMC, and that the 
metal-rich stars appear to have relatively long periods for metal-rich RRab stars (i.e., 
the subgroup of RRab stars with shorter periods at a given c^). 
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4. Physical Characteristics based on Fourier Parameters 

Several different physical characteristics were established using equations derived 
by Jurcsik (1998, hereafter J98), Kovacs & Walker (2001), and Sandage (2004, 2006). 
Since these equations were originally established for V data, and we have only Kp data, 
it was necessary first to establish Kp versions of the equations. The specific equations 
used to derive the physical quantities, and the detailed tables given in Nemec et al. 
(2011) are not repeated here, but the results of their application are summarized in 
several graphs and discussed below. 

Figure 6 shows four different metallicity diagnostic diagrams. The upper-left panel 
is the classical period-amplitude diagram, often called the 'Bailey P-A diagram' but 
coined the 'Oosterhoff-Arp-Preston P-A diagram' by Sandage (2004, hereafter S04). 
At a given total-amplitude lower [Fe/H] stars tend to have longer periods. Spectro- 
scopically calibrated diagonal lines, such as that shown in the lower right panel for the 
Messier 3 RRab stars = 3.124 + 5.1281ogP), can be constructed, from which offsets 
give [Fe/H] values. Since the pattern of the stars seen in all four panels is very similar 
all four variables can, in principle, be used to derive [Fe/H] values. Four different es- 
timates of [Fe/H] were calculated by Nemec et al. (2011), based on the ^ relation of 
J98 (her equation 1), the 0^ relation of S04 (his equation 3), the A tot relation of S04 
(his equation 6), and the 'risetime' relation of S04 (his equation 7). The 03i relations 
seem to be the most reliable and the [Fe/H] values based on the J98 formula, after 
conversion to the Zinn-West scale, are given in the last column of Table 2. Recently 
acquired CFHT 3.6-m ESPaDOnS spectra confirm that KIC 6100702 is indeed metal- 
rich, at [Fe/H] = -0.18±0.06 dex, while AW Dra is of intermediate metallicity, with 
[Fe/H] = -1.33±0.08 (Nemec, Ripepi, Chadid et al. 2011, in preparation). 

Figure 7 summarizes the derived de-reddened colour information in the form of 
correlation diagrams. In all four panels the four more metal rich stars stand apart from 
the more metal-poor stars. For this reason they have been plotted with open squares. 
The linear correlations show that the reddest metal-poor stars (those nearest the red 
edge of the Instability Strip) are more metal poor, have longer periods and smaller total 
amplitudes, and are more luminous than the bluest metal-poor stars. The diagrams 
also show that at a given colour the four metal-rich stars (plotted with open squares) 
have shorter periods, smaller amplitudes and are less luminous than the metal-poor 
stars. 

Not represented in these diagrams but computed and given in the full paper are 
absolute magnitudes (V-passband), surface gravities, mean effective temperatures, dis- 
tances, and luminosities and masses. Both pulsational and evolutionary L and M. val- 
ues were computed. The L(puls) were calculated using equations 16 and 17 of J98, and 
in both cases the lower metallicity stars have the higher L(puls). The A^(puls) were 
calculated using equations 14 and 22 of J98. The average mass for the four metal-rich 
stars is ~0.50 Ai Q compared with the average mass for the more metal-poor stars 
~0.60 M.q. Three different L(evol) were calculated, using equations 8, 10 and 12 of 
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Figure 6. — Period- amplitude and three other metallicity diagnostic diagrams (all with log P along 
the abscissa) for the Kepler non-Blazhko RRab stars. AW Dra, FN Lyr and NR Lyr are plotted twice 
in each panel: the Kepler points are shown as black squares, the T^-photometry points are shown as 
open diamonds, and the two are connected by vertical lines. The consistency of the Kp-V offsets seen 
here and elsewhere justifies using the J98 and Sandage relations to calculate physical characteristics. 
In the lower-right panel the four metal rich stars have been plotted as open-squares, and the dashed 
diagonal line is the observed ^-relation for the M3 RR Lyrae stars (from Cacciari et al. 2005). 



Sandage (2006). For the metal-poor stars all three L(evol) are systematically larger 
than the L(puls) values. For the four metal-rich stars the agreement is better but 
there is a wide range of L(evol) owing to the linear or non-linear [Fe/H] dependencies. 
Regardless of which formula was used the most luminous stars have the lowest metal- 
licities and there is internal consistence. It is not clear whether the L(puls) or L(evol) 
are correct. As was the case for the luminosities, the A^(evol) are all larger than the 
corresponding M. (puis) values. 

Figure 8 compares the results with the locations of model zero-age horizontal 
branches and evolutionary tracks. The top two panels show HR diagrams with log T e s 
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Figure 7. — Four diagrams for the non-Blazhko RR Lyr stars, all with mean dereddened colour 
(B — V)o (average of the two values given in column 6 of Table 5) along the abscissa: metallicity vs. 
colour (Top left), period vs. colour (Top right), total amplitude vs. colour (Bottom left), and absolute 
magnitude vs. colour (Bottom right) . In each panel the line (and its equation) is from a least squares 
fit to the points for the 15 metal-poor stars (i.e., those with [Fe/H]< —1.0 dex). 



as abscissa and logL as ordinate, and two sets of ZAHB loci from Dorman (1992) 
computed for different masses along the ZAHB. The more luminous horizontal branch 
assumes [Fe/H]= —1.78 dex, and the less luminous branch [Fe/H]= —0.47 dex. The 
numbers next to the symbols are the masses for the individual models, which are seen to 
be higher for the low-[Fe/H] tracks than for the high-[Fe/H] tracks. For both assumed 
metallicities oxygen enhanced and non-enhanced ZAHBs have been plotted - the effect 
of increasing the oxygen to iron ratio from [O/Fe]=0 to 0.66 for the low-metallicity 
ZAHBs is to lower the luminosity and reduce the mass at a given temperature. For the 
high-metallicity tracks an oxygen enhancement from [O/Fe]=0 to 0.23 has little effect 
on the derived L or M.. Also plotted in the low-[Fe/H] case are the evolutionary paths 
away from the ZAHB for two masses, 0.66 and 0.68 M. Q . In both panels the non- 
Blazhko RR Lyr stars with low metallicities are represented by large black squares, the 
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Figure 8. — (Top panels) HR-diagrams showing the locations of the Kepler non-Blazhko RR Lyr 
stars compared with two sets of horizontal branch models (low- and high-metallicity) from Dorman 
(1992). In both graphs the four metal-rich RR Lyr stars have been plotted with open squares and 
the metal-poor stars with filled squares; the black numbers next to the ZAHB symbols and the two 
evolutionary tracks are the assumed masses for the individual stellar evolution models. For the non- 
Blazhko RR Lyr stars in the top left panel the L(puls) were calculated with eq. 17 of J98 and labelled 
(in red) with the .M(puls) calculated with eq. 22 of J98; in the top right panel the L(evol) were 
calculated with eq. 12 of S06 and labelled (in red) with the .M(evol) calculated with eq. 15 of S06. 
Both panels (left and right) show the blue and red edges for the instability strip computed using 
the Warsaw pulsation code; the fundamental mode edges are plotted as solid lines, and first-overtone 
mode edges as dashed lines. (Bottom panels) Graphs comparing the observed pulsation periods and 
the pulsation periods derived with the Warsaw pulsation code. The assumed masses and luminosities 
are the values given in the top panels. For each star two points are plotted, computed assuming the 
T e g from eq. 11 of KW01 (blue dots) and from eq. 18 of S06 (red circles). The best agreement is 
seen in the left panel, with the KW01 temperatures favoured over the S06 values. In both panels the 
largest differences are seen for the longest period stars, with no metallicity dependence. 

four high-metallicity stars are plotted with open squares, and the T e s are the average 
of the Kovacs & Walker (2001, hereafter KW01) and Sandage (2006, hereafter S06) 
values. 

In the top left panel of Fig. 8 the luminosities and masses (labelled in red) of 
the Kepler non-Blazhko RR Lyr stars were calculated with eqs. 17 and 22 of J98 
and thus are based on pulsation theory. The L(puls) and .M(puls) are seen to be 
systematically smaller than values derived from the ZAHB tracks (for the appropriate 
metal abundance). The reddest non-Blazhko RR Lyr stars lie close to the fundamental 
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mode red-edge and have the smallest amplitudes. This graph also shows blue and red 
edges of the instability strip for the fundamental mode (red and blue solid lines) and 
first-overtone mode (red and blue dashed lines). The edges were calculated with the 
Warsaw pulsation code assuming a mass of 0.65 M Q . The Kepler non-Blazhko stars all 
lie in the fundamental mode region of the variability strip, and the smallest amplitude 
RR Lyr stars (the four metal-rich stars, KIC 7030715 and NR Lyr) have locations near 
the fundamental red edge (FRE) of the instability strip. As expected, all the stars near 
the FRE have low R 3i values. 

In the top right panel of Fig. 8 the luminosities were calculated with eq. 12 of S06, 
and the masses (labelled in red) with eq. 15 of S06; thus they are evolutionary L and 
M. values. In this case there is very good agreement with the stellar evolution models, 
as one expects since they are based on stellar evolution models. Enhancing the oxygen 
to iron ratio by the plotted amounts makes little difference. 

It is unclear which are correct, the L(puls) and .M(puls), or the L(evol) and 
.M(evol). The mass and luminosity discrepancies go in the same direction as seen 
for Cepheids. Pietrzynski et al. (2010) recently derived a dynamical mass .M(dynam) 
for a classical Cepheid in a well detached, double-lined eclipsing binary in the LMC. 
The mass they derive is very accurate and favours A^(puls). The reason for the dis- 
crepancies may be the same, as suggested by Pietrzynski et al. - not enough mass loss 
has been taken into account in the evolution models. 

Finally, the Nemec et al. (2011) paper concludes with a section on hydrodynamic 
models computed with the Smolec and Moskalik 'Warsaw convective pulsation pro- 
gram'. Perhaps most relevant for the discussion above is the conclusion that in the 
P-A diagram when [Fe/H] is varied from —0.2 to —1.9 dex there is very little motion 
at a given amplitude. Instead, motion to longer periods at a given amplitude occurs 
for stars of higher L for a given mass, or for stars of lower mass for a given luminosity. 
This suggests that the main factors causing period shifts in the P-A diagram seem to 
be the mass and luminosity, the latter being caused either by higher ZAHB mass or by 
post-ZAHB evolution. 

5. Conclusions 

The main conclusions reached by Nemec et al. (2011) concerning the Kepler non- 
Blazhko RR Lyrae stars are: (1) Fourier-based [Fe/H] values have been derived for the 
stars (using four different correlations), and four of the 19 stars (V782 Cyg, V784 Cyg, 
KIC 6100702 and V2470 Cyg) appear to be much more metal-rich than the others - 
recent CFHT spectra confirm that KIC 6100702 is metal-rich with [Fe/H]= —0.18 ± 
0.06 dex, and AW Dra is of intermediate metallicity with [Fe/H]=— 1.33 ± 0.08) dex. 
(2) The star KIC 7021124 is found to be doubly-periodic (P 2 /P = 0.59305) with 
characteristics similar to V350 Lyr (see Benko et al. 2010); (3) All the stars show 
remarkably constant light curves (a < 1 mmag) over the 420-day interval of the Q0- 
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Q5 observations; (4) FN Lyr and AW Dra are found to have increasing periods; (5) 
the ASAS-North photometry of nine of the brightest stars shows that their colours 
are bluest when brightest, as expected for RR Lyrae stars; (6) Fourier-based physical 
characteristics for all the stars have been derived; (7) Distances range from 3 to 23 Kpc; 
(8) Pulsational and evolutionary masses and luminosities are calculated for the stars, 
the latter tending to be higher in both cases. When the observed periods are compared 
with periods computed with the Warsaw non-linear convective pulsation code better 
agreement is achieved assuming pulsational L and M values rather than the (higher) 
evolutionary L and M values. (9) Finally, the Warsaw models show that varying [Fe/H] 
has a relatively small effect on the P-A relation and that the main factors causing the 
period shifts must be luminosity and mass. 
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